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Grand unification of AGN activity in the ACDM cosmology 
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ABSTRACT 

We track the coevolution of supermassive black holes (SMBHs) and their host galaxies 
through cosmic time. The calculation is embedded in the GALFORM semi-analytic model 
which simulates the formation and evolution of galaxies in a cold dark matter (CDM) uni- 
verse. The black hole (BH) and galaxy formation models are coupled: during the evolution 
of the host galaxy, hot and cold gas are added to the SMBH by flows triggered by halo gas 
cooling, disc instabilities and galaxy mergers. This builds up the mass and spin of the BH, 
and the resulting accretion power regulates gas cooling and subsequent star formation. The 
accretion flow is assumed to form a geometrically thin cool disc when the accretion rate ex- 
ceeds O.OlAfEdd, and a geometrically thick, radiatively inefficient hot flow when the accretion 
rate falls below this value. The resulting quasar optical luminosity function matches observa- 
tions well, and the mass of the SMBH correlates with the mass of the galaxy bulge as in the 
observed Mbh — Mbuige relation. The BH spin distribution depends strongly on whether we 
assume that the gas in any given accretion episode remains in the same plane or it fragments 
into multiple, randomly aligned accretion episodes due to its self-gravity. We refer to these 
cases as the "prolonged" and "chaotic" accretion modes respectively. In the chaotic accre- 
tion model there is a clear correlation of spin with SMBH mass (and hence host galaxy bulge 
mass). Massive BHs (M > 5 x 10 s M©) are hosted by giant elliptical galaxies and are rapidly 
spinning, while lower mass BHs are hosted in spiral galaxies and have much lower spin. Us- 
ing the Blandford-Znajek mechanism for jet production to calculate the jet power, our model 
reproduces the radio loudness of radio galaxies, LINERS and Seyferts, suggesting that the 
jet properties of active galaxy nuclei (AGN) are a natural consequence of both the accretion 
rate onto and the spin of the central SMBH. This is the first confirmation that a CDM galaxy 
formation model can reproduce the observed radio phenomenology of AGN. 
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1 INTRODUCTION 

Active galaxies can be classified according to the importance of the 
radio emission from their nucleus. Objects showing strong emis- 
sion at radio wavelengths belong to the "radio-loud" class, whereas 
those with negligible emission belong to the "radio-quiet" class. 
Radio-loud AGN are associated with large-scale radio-emitting jets 
while radio-quiet AGN show very little or negligible jet activity. 
Sikora et al. (2007) found that AGN have a bimodal distribution 
on the radio-optical plane, with radio-loud objects being about I0' 5 
times brighter in radio than radio-quiet objects (see also Keller- 
mann et al. 1989; Xu et al. 1999). The origin of this dichotomy 
remains unknown. However, it has been proposed in many studies 
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that a first step towards explaining why some AGN launch promi- 
nent jets while others do not is to understand the nature of the cen- 
tral engine, the accreting BH. 

AGN jets are believed to extract rotational energy from the 
BH and accretion disc through magnetic fields. Analytical stud- 
ies have contributed to an understanding of the nature of the jets 
(Blandford & Znajek 1977; Macdonald & Thorne 1982; Blandford 
& Payne 1982; Begelman, Blandford & Rees 1984), but a break- 
through has come from magneto-hydrodynamical (MHD) simu- 
lations of the accretion flow. These self-consistently calculate the 
turbulent magnetic field dynamo which is the physical origin of 
the stresses which transport angular momentum outwards, allow- 
ing material to accrete onto the BH (Balbus & Hawley 1998). Em- 
bedding such calculations in a fully general relativistic framework 
produces relativistic jets from the accretion flow, and the collima- 



© 0000 RAS 



2 Fanidakis et al. 



tion and jet power depend on the BH spin (McKinney & Gammie 
2004; De Villiers et al. 2005; Hawley & Krolik 2006). 

Thus, it now seems clear that the jet power depends on the 
vertical (poloidal) magnetic field strength close to the BH, but the 
simulations are still not yet at the level where they can predict this 
ab. initio (e.g. Beckwith et al. 2008), so analytic models are still 
required. The most popular of these, the Blandford-Znajek (here- 
after BZ) mechanism, uses the magnetic field as a means to tap 
the spin of the BH. This gives a strong spin dependence on the jet 
power. Indeed, many authors have proposed spin as the physical pa- 
rameter that determines the radio loudness of an AGN (see Wilson 
& Colbert 1995; Hughes & Blandford 2003). The spin paradigm, 
as this idea is called, offers a plausible theoretical explanation for 
the wide range of jet luminosities in AGN and could be the basis 
for understanding the observed dichotomy between radio-loud and 
radio-quiet AGN. 

Alternatively, the radio-loud, radio-quiet switch may just be 
determined by the physical state of the accretion flow. Stellar-mass 
BH binary systems in our galaxy show a clear spectral transition at 
about 1% of the Eddington luminosity, L ~ 0.01Z/Edd, from a hot, 
optically thin accretion flow at low luminosities (e.g. an advection- 
dominated accretion flow or ADAF; Narayan & Yi 1994) to a cool, 
geometrically thin disc (Shakura & Sunyaev 1973) at higher lumi- 
nosities (Esin et al. 1997, see e.g. the review by Done, Gierlinski 
& Kubota 2007). The collapse of the radio jet observed across this 
transition clearly relates to the large drop in pressure (and hence 
scale height, H) between the hot and cool flow (e.g. Fender, Bel- 
loni & Gallo 2004). Even without BH spin, this produces a clear 
dichotomy of radio properties (e.g. Jester 2005). 

It seems most plausible that both these mechanisms, along 
with the mass accretion rate, affect jet power. The BZ mechanism 
has a "hidden" dependence on mass accretion rate and the scale 
height of the flow because the magnetic field strength close to the 
BH saturates, due to the dynamo, into rough equipartition with the 
pressure in the flow. This depends on the mass accretion rate for ei- 
ther the disc or the hot flow, but the pressure in the hot flow is much 
larger than that in the cool disc so the field strength and hence jet 
power abruptly drop at this transition (Meier 2001; 2002) 

Sikora et al. (2007) used these ideas to interpret observations 
of radio galaxies and concluded that the data could be explained if 
the jet depends on mass accretion rate, accretion model and spin. 
However, they also required a mass-spin correlation, such that the 
most massive BHs, which, according to the Mbh — Mb u i ge relation- 
ship (Magorrian et al. 1998; McLure & Dunlop 2002; Marconi & 
Hunt 2003; Haring & Rix 2004), reside in massive elliptical galax- 
ies, have higher spin than lower mass BHs, which reside in spi- 
rals. They speculated that this could occur through the hierarchical 
growth of structure, where the major mergers which give rise to gi- 
ant elliptical galaxies trigger large amounts of gas accretion with 
a given angular momentum direction onto the BH, spinning it up 
to the maximal value. By contrast, spiral galaxies have not expe- 
rienced a recent major merger, and grow mainly through smooth 
accretion and multiple minor mergers with random angular mo- 
mentum, resulting in low spin and hence weak jet luminosities. 

Some of these ideas have recently been explored using simpli- 
fied models of galaxy formation by Berti & Volonteri (2008) and 
Lagos et al. (2009). In this paper, we calculate the growth of BHs 
by accretion and mergers, their acquisition of spin and their accre- 
tion rates within a specific theory of galaxy formation in a CDM 
model which has been extensively tested against a large range of 
observations of the galaxy population. Not only is our calculation 
of the joint evolution of BHs and galaxies self-consistent, but we 



also incorporate several accretion and jet-launching models from 
the literature which allows us to perform a quantitative comparison 
with observations of AGN. 

The paper is organised as follows. In Section 2, we explain 
how we calculate the growth of the mass and spin of SMBHs 
through the accretion of hot gas from the halo (radio mode) and 
cold gas from the galactic disc and mergers (quasar mode). Most 
BHs grow primarily by gas accretion except the most massive ones 
which form late and increase their mass substantially by merging 
with other SMBHs. A BH is spun up to approximately the maximal 
value when it accretes its own mass from a flow at constant angular 
momentum. 

In Section 3, we present the resulting distributions of BH 
mass, gas accretion rate and spin. We consider two possible modes 
of accretion. In the "prolonged accretion" case, the gas is accreted 
in a single episode. Even minor mergers trigger gas flows onto the 
nucleus that often deposit a mass greater than the mass of the recipi- 
ent BH. Thus, most BHs are typically spun up to the maximal value. 
This model predicts radio properties that do not reproduce the ob- 
served luminosity function. We then consider a "chaotic accretion" 
case in which the accretion episodes are limited by the self-gravity 
of the disc. Gas flows onto the nucleus give rise to a series of ac- 
cretion episodes each typically augmenting the BH mass by only a 
small factor. Successive accretion events are uncorrelated, resulting 
in low spins for the relatively low mass BHs which grow primarily 
by accretion. On the other hand, the most massive holes which build 
up substantial mass through BH-BH mergers are spun up to high 
(but not maximal) spin values. In Section 4, we show that adopting 
this accretion model, the optical luminosity from disc-accreting ob- 
jects matches the quasar luminosity function reasonably well. 

In Section 5, we incorporate an explicit BZ model for the jet 
power (Meier 2002) and use this to predict the radio luminosity 
function which, in the chaotic accretion case, agrees well with ob- 
servations. On the basis of our derived optical and radio luminos- 
ity distributions, we propose, in Section 6, a "grand unification of 
AGN activity", in which the accretion flow and jet luminosity are 
related to AGN type through their mass, spin and mass accretion 
rate. Fundamentally, our calculation shows, for the first time, how 
the coeval growth of BHs and their host galaxies result in optical 
and radio properties that explain the AGN activity seen in the local 
Universe. 



2 EVOLUTION OF SPIN IN HIERARCHICAL GALAXY 
FORMATION MODELS 

In this Section we describe the cosmological processes that initi- 
ate and regulate the growth of SMBHs and outline the modelling 
of the spin evolution due to gas accretion and mergers in hierar- 
chical cosmologies. Our basic modelling tool for our predictions 
is the GALFORM code (Cole et al. 2000), and we use an update of 
the Bower et al. (2006) version for modelling the formation and 
evolution of galaxies in the ACDM cosmology. The changes rela- 
tive to the Bower et al. model are the following. Firstly, the frac- 
tion, esmbh, of the Eddington luminosity of an accreting SMBH 
that is available for heating the halo during an episode of AGN 
feedback is set to 0.01 (Bower et al. 2006 use, esmbh = 0.o43). 

1 Note that due to an error in Bower et al. (2006), cooling luminosities were 
overestimated by a factor of 4tt. Thus, while the paper quotes the efficiency 
parameter esMBH as 0.5, this should have been 0.5/47T = 0.04. With this 
correction, the rest of the parameters and results are unchanged. 
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Secondly, in starbursts triggered by galaxy mergers or disc insta- 
bilities, we assume that the fraction, Fsmbh, of the cold gas turned 
into stars that is accreted onto the BH is 0.01 (Bower et al. 2006 
use Fsmbh = 0.017). These changes are introduced to improve 
the modelling in this paper and do not affect the fundamental pre- 
dictions of the Bower et al. model. 

BHs are assumed to evolve in mass in accordance with the 
model developed in Malbon et al. (2007) and Bower et al. (2006). 
We study the cosmological spin evolution of SMBH seeds of a to- 
tal of ~ 4.2 x 10 6 galaxies identified in the Millennium N-body 
simulation (co-moving volume of 1.25 x 10 8 fr~ 3 Mpc 3 ) from red- 
shift 127 to redshift zero (Springel et al. 2005). Here h is defined 
by Ho = h x 100 km s _1 Mpc _1 , where Ho is the Hubble con- 
stant at redshift zero. The cosmology adopted in the simulations is 
h = 0.72, Q m = 0.25, fl h = 0.045, fl A = 0.75 and cr 8 = O.fl 

2.1 The growth of SMBHs 

The evolution of BH mass fits naturally in the CDM model 
of galaxy formation (Kauffmann & Haehnelt 2000; Malbon et 
al. 2007), where structures grow hierarchically. Small structures 
form first, then evolve through mergers into large ones. In the event 
of a galaxy merger, the less massive galaxy (satellite) sinks into the 
gravitational potential of the massive central galaxy as a result of 
dynamical friction (Binney & Tremaine 1987). In our model we as- 
sume that if the mass of the satellite galaxy is comparable to that of 
the central galaxy, the merger disrupts the galaxies and resulting in 
the formation of an elliptical galaxy (major merger). This event is 
accompanied by a burst of star formation as the available cold gas 
from both progenitors is transferred to the centre and transformed 
into stars. Some of this cold gas reservoir feeds the central SMBH. 
An additional process of cold gas accretion that contributes to the 
SMBH mass is the collapse of galaxy discs triggered by dynami- 
cal instabilities (Efstathiou, Lake & Negroponte 1982). When the 
self-gravity of the galactic disc becomes sufficiently large, the disc 
equilibrium is disrupted resulting in the formation of a bar which 
enables gas to be transferred towards the centre as the disc. A frac- 
tion of that gas is directly fed into the BH through an accretion 
disk and powers an AGN (Lynden-Bell 1969), while the rest un- 
dergoes star formation. For highly efficient accretion activity, the 
disc/BH system dominates the energetics of the nucleus, and the 
galaxy becomes visible as a quasar. SMBHs at the centres of galax- 
ies observed in the local Universe are regarded as the remnants of 
quasar phases at earlier epochs. 

The consequence of a minor merger, namely, the accretion of 
satellites of low mass compared to the central galaxy, is usually 
less dramatic. During a minor merger in our model, cold gas and 
stellar content of the satellites are added to the central galaxy, and 
a merger-driven starburst may supply fresh gas to the central BH. 

SMBHs acquire part of their mass merging with other 
SMBHs. The formation of a BH-BH binary and the subsequent co- 
alescence of the BHs is a natural evolutionary stage for a SMBH if 
the host galaxy experiences a merger. As the galaxies merge, we as- 
sume that the SMBHs sink to the centre by dynamical friction from 
distant stars or by viscous effects from the surrounding gas. The 
transition to a bound binary state after the galaxy merger is an open 

2 f2 m , Hj, and SI a express the present density of the baryonic, total matter 
and dark energy components of the Universe relative to the critical density 
(Pcrit = 3H 2 /8nG). erg measures the rms mass fluctuations in spheres of 
radius 8 ri -1 Mpc linearly extrapolated to the present day. 



issue (Milosavljevic & Merritt 2001). However, it is believed that 
once the separation of the two BHs becomes small enough, gravita- 
tional radiation carries away the remaining angular momentum of 
the binary. The removal of energy from a SMBH binary due to grav- 
itational wave emission leads to a gradual shrinkage of the relative 
separation of the two members. There is a point where the eccen- 
tricity reaches zero and the orbit circularises. At that time, the two 
BHs are very close to each other, and gravitational waves are emit- 
ted copiously. The radiated energy is so large, that two SMBHs, 
which are a few AU apart, lose all their potential energy within a 
couple of minutes and inevitably coalesce. After the coalescence is 
complete the binary enters the ringdown phase, where the merged 
members settle into a quiescent remnant hole. 

A third channel for SMBH growth in our model is provided 
by diffuse gas in dark matter haloes undergoing quasi-hydrostatic 
cooling. When a massive halo collapses gas is shock heated out to a 
radius comparable to the virial radius of the dark matter halo. These 
haloes have a cooling time that is longer than the free-fall time of 
the gas and, thus, the gas settles into a quasi-static atmosphere sur- 
rounding the galaxy rather than simply falling towards the centre 
(White & Frenk 1991). This atmosphere - the "hot halo" regime 
- is in pressure supported hydrostatic equilibrium and extends be- 
yond the virial radius of the dark matter halo. The galaxy is then 
supplied with cold gas by cooling flows at the disc centre, which 
also feed the central SMBH. 

The formation of massive hot haloes would lead to the growth 
of very massive galaxies unless a heating mechanism regulates the 
cooling flow. Bower et al. (2006) and Croton et al. (2006) invoke 
energy injection from the central SMBH during the so-called "radio 
mode" feedback (see also de Lucia et al. 2006; Cattaneo et al. 2007; 
Lagos et al. 2008). The radio mode is assumed to occur during the 
quiescent accretion of gas from the hydrostatically supported hot 
halo onto the SMBH, during which energy from the BH accretion 
is injected directly into the hot halo suppressing the cooling flow. 
In the Bower et al. (2006) model, the cooling flow stops when the 
power from the SMBH is sufficient to offset the rate at which en- 
ergy is being radiated away. 

Bower et al. (2006) adopt a BH growth model in which during 
a disc instability or galaxy merger, the BH accretes a fixed frac- 
tion of the gas that turns into stars in the burst taking into account 
processes such as feedback and recycling in the galaxy (Malbon 
et al. 2007). The amount of gas deposited onto the BH is set by 
an efficiency factor, which determines the fraction of the available 
gas reservoir for star formation that is accreted by the hole. The 
value of that parameter is chosen to fit the normalisation of the local 
A/bh — A/buige relation. Hereafter, we refer to the accretion of cold 
gas triggered by disc instabilities and galaxy mergers as "quasar 
mode" accretion and to the accretion from quasi-hydrostatic haloes 
as "radio mode" accretion, following the terminology introduced 
by Croton et al. (2006). 

2.2 Astrophysical processes affecting BH spin evolution 

SMBHs are expected to possess angular momentum Jbh = 
\a\GM^ h /c, where a is the spin parameter, ^ |o| ^ 1. The spin 
has a significant impact in the close vicinity of the BH. For exam- 
ple, it determines the efficiency for converting matter into radiation 
in an accretion disc (Novikov & Thorne 1973) and it is believed 
to influence the formation and direction of the radio jets in AGN 
(Blandford & Znajek 1977; Macdonald & Thorne 1982; Begelman, 
Blandford & Rees 1984). In addition, it is of special interest in the 
modelling of gravitational waves from BH-BH binaries (Merritt et 
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al. 2005; Baker et al. 2006a, 2007; Campanelli et al. 2007; Buo- 
nanno et al. 2007). 

The evolution of spin is closely related to the channels that 
contribute to the growth of the BH. Each mechanism for BH growth 
is associated with different spin evolution. For example, accretion 
of gas that co-rotates with the BH should spin up the hole (Bardeen 
1970), whereas, the merger of two equal-mass non-spinning BHs 
results in a final spin of 0.69 (Baker et al. 2006b; Berti et al. 2007, 
Hinder et al. 2008). We explain below how these mechanisms that 
influence the spin of a BH are included in our model. 

2.2.7 Gas accretion 

After the SMBH seed forms at the centre of a galaxy, accretion 
usually initiates the growth era. We assume that an accretion disc is 
formed in the equatorial plane of the hole. As proposed by Lynden- 
Bell (1969), gas parcels in the disc gradually lose angular momen- 
tum due to viscous torques exerted by magnetic fields and drift ra- 
dially inwards until they reach the inner edge of the accretion disc. 
The inner edge of the disc is usually taken as the location of the 
last stable orbit (LSO) around the BH. The LSO is a function of 
the hole's angular momentum and can be written as (Bardeen et 
al. 1972): 

n 8 o/-Rg = {3 + Z 2 ±[{3-Z 1 ){3 + Z 1 +2Z 2 )] 1/2 }, (1) 

is defined 

and Zi, Z 2 are defined in terms of the spin, a, as 



'''Is 

where the gravitational radius, R g 
Schwarzschild radius of the BH, i?schw 



half of the 

2GM hh /c 2 , 



Zi 

Z 2 



1 + 



(l-a 2 ) 1/3 [(l + a) 1 / 3 + (l-a) 1 / 3 ] 



(3a 2 + Zl) 1/2 . 



(2) 



Note that when the spin parameter has a negative sign the BH is 
counter-rotating with respect to the orbit of a particle around it. 
Then, for counter-rotating orbits (— 1 a < 0) 9 ^ n so < 6 and, 
for co-rotating orbits (0 < a ^ 1), 6 < n so ^ 1. 

An important property is the binding energy of the gas at the 
LSO, defined as the difference between the rest energy of a gas par- 
cel at infinity and its energy at the LSO as measured by an observer 
at infinity. If e expresses the energy per unit rest mass, we can de- 
fine the binding energy as 1 — e/c 2 . This provides a simple relation 
for the accretion efficiency, 



ei s 



, 2 

d/C 



1 



(3) 



_ 2 1 

3 f i so ' 

(Novikov & Thorne 1973) which corresponds to the fraction of the 
energy released by matter spiralling in towards the BH through a 
succession of almost circular orbits. Eq. l[3) shows that accretion of 
matter onto slowly rotating BHs has moderate efficiency. For non- 
rotating BHs (a = 0, n so = 6) this is just 1 - or 5.7%. 
However, for co-rotating matter around rapidly rotating BHs, the 
efficiency increases significantly reaching 1 — l/\/3 or 42.3% as 
a — > 1. This sets an upper limit to the efficiency of the accretion 
onto BHs, a process obviously much more efficient than thermonu- 
clear burning. 

Once the gas reaches the LSO, we assume that it falls directly 
into the BH. In this way, the accretion carries into the BH the energy 
per unit mass, e, and angular momentum per unit mass, /, that the 
gas has at the LSO. Thus, accretion of a rest mass dMo leads to a 
change in the hole's total mass, Mbh, and angular momentum, Jbh, 
equal to, 
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Figure 1. The final BH spin after the accretion of gas in a counter- and 
co-rotating configuration. The dotted lines indicate the final mass needed to 
spin up a non-rotating and maximally counter-rotating BH to a? = 0.998. 



The change in the hole's spin induced by the accretion of dMo is 
governed by the differential equation 



da 



Ik 



din M b 



Mbh G eiso 



2a. 



(5) 



This was integrated by Bardeen (1970) using the explicit expres- 
sions for e and I, resulting in the following solution, 



f _ 1 .1/2 M bh 
a ML 



3 a 



Mbh 
ML 



- 2 



1/2- 



(6) 



when M^/ M hh 



< f!/ 2 . Here, M* 



and are the final mass 



dMbh = (ei S o/c )dM , dJ b h = /i so dM . 



(4) 



and spin parameter of the BH. If ML/M^h > r^J 2 , then the final 
spin is always equal to unity. The expression in Eq. @ governs 
the evolution of a during accretion from an initial state of a co- 
rotating or counter-rotating disc. According to this, a non-rotating 
BH will be spun up to a maximum rotation (a = 1, fi so = 1) 
after increasing its mass by 1.44Mbh (see Fig.[TJ. For a maximally 
rotating BH in a counter-rotating accretion disc (a = — 1, ri 8 o = 
9), an increase in mass of 2Mbh is required in order for the BH 
to be spun down to a = and subsequently spun up to maximum 
rotation. As implied by Eq. l[6j, further accretion onto the BH keeps 
a equal to unity. 

Bardeen's calculations are limited to the process of angu- 
lar momentum transport between the accreted matter and the BH, 
without taking into account any other processes. In fact, Thorne 
(1974) argued that the accretion disc radiates and that some of this 
radiation will be accreted by the hole. Capture of photons with 
angular momentum opposite to that of the BH will then produce 
a counteracting torque that prevents spin up beyond the limiting 
value of a = 0.998 (fi so = 1.23, e ~ 0.32). In addition, axial 
relativistic outflows of matter in the form of jets may have strong 
implications for the upper spin limit, since these outflows are ac- 
celerated by magnetic fields that are powered by the extraction of 
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Accretion Disk 




Figure 2. Schematic illustration of a warped accretion disc. Jbh is the an- 
gular momentum of the BH (| Jbhl = \ a \GM^ h /c), is the angular 
momentum of the disc given by Eq. \2Q\ and Jtot represents the total an- 
gular momentum of the system, ,/bh + Jd ■ 



the rotational energy of the BH (Hawley & Krolik 2006; Benson & 
Babul 2009). 



2.2.2 Gas accretion through a misaligned disc 

In the general case where the accretion disc does not lie on the 
equatorial plane of the BH but has a random orientation rela- 
tive to the angular momentum of the BH, the evolution of the 
spin is a complicated process. Following the discussion in King 
et al. (2005), we assume that a thin disc is inclined at some random 
angle 8 relative to the orientation of the hole's angular momentum 
vector, Jbh- We denote the angular momentum of the disc as Jd 
(see King et al. 2005 and Volonteri et al. 2007 for a discussion of 
the nature of Jd), and define the total angular momentum vector 
J tot as 



= Jbh + Jd 



(7) 



The angle 8 between Jbh and Jd is defined such that ^ 8 ^ n; 
the values 8 = and 8 — tt correspond to full alignment and 
anti-alignment respectively. J to t is a constant vector i.e. has fixed 
orientation as well as magnitude for a given accretion event. Its 
magnitude is given by 



Jtot — Jbh 



2JbhJd cosf?. 



(8) 



We further introduce the angle 8 t as the angle between Jbh and 

J tot- 

When the vectors Jbh and Jd are misaligned, the tilted or- 
bits of the gas parcels in the accretion disc experience a torque due 
to the Lense-Thirring effect, which causes the plane of the accre- 
tion disc to precess about the rotational axis of the BH (Lense & 
Thirring 1918; Wilkins 1972; Bardeen & Petterson, 1975; Scheuer 
& Feiler, 1996). If viscosity is strong enough this can force the in- 
ner parts of the disc to rotate into the equatorial plane of the hole 
resulting in a "warped disc" (see Fig.[2]l. 

The Lense-Thirring torque can be expressed as 



dL 

dt 



rip x l, 



(9) 



where L is the angular momentum per unit area of the disc and tt p 
is the precession rate and is given by 



S7, 



2G J bh 

c 2 R? 



where R is the distance from the BH (Pringle 1992). The precession 
timescale is therefore defined as, 



2tt 



(11) 



PIUC " n p (R)' 

which is proportional to R 3 and thus is much shorter closer to 
the BH. Other timescales relevant to this problem are the viscous 
timescales of accretion, t vi , and warp propagation, t„ 2 , 

= ^ry (12) 

Here, v\ and v-i are the kinematic viscosities acting on velocity gra- 
dients parallel and normal to the plane of the disc respectively. The 
balance between the timescales t plcc and t V2 determines whether 
the Lense - Thirring torque is able to align the inner disc with the 
spin axis. The condition for alignment is t prcc < t„ 2 (Natarajan & 
Armitage 1999). In other words, the disc will be aligned with the 
spin at radii where the precession timescale is much shorter than the 
timescale of radial diffusion of the warp. The characteristic radius, 
fiwarp, of the aligned part of the disc follows from the condition 

tprec — ' 

The evolution of Jbh and Jd during the precession is deter- 
mined by the conditions (King et al. 2005): 



d J 2 -0 



-ijj<0 and -icosft^O. (13) 
dt at 



As can be inferred by the first two conditions, the magnitude of 
Jbh remains constant while the angle 8 t decreases with time, which 
implies that Jbh always aligns with J to t ■ In contrast, Jd decreases 
in magnitude as Jbh aligns, which is to be expected, since the total 
angular momentum has to remain constant. 

The end result of the Lense - Thirring precession is a BH 
which is aligned or anti-aligned with the surrounding accretion 
disc. The expression for the magnitude of J to t allows us to exam- 
ine the final configuration of the system. Obviously, if J^ h > J t 2 ot 
anti-alignment occurs, which requires 



cos 8 < 



Jd 

2Jbh 



(14) 



Hence, a BH with 8 > n/2 and 2 Jbh > Jd, eventually anti-aligns 
with the accretion disc. This means that if Jd > 2 Jbh then align- 
ment always happens as the disc angular momentum completely 
overwhelms that of the BH. 

Further investigation of the warped discs requires knowledge 
of the accretion flow properties. We use the thin disc solution 
of Shakura & Sunyaev (1973) for our analysis. In the standard 
Shakura-Sunyaev disc model, the analytic expression for the warp 
radius depends on the values of Jbh, AJbh, and the accretion rate M 
of the BH. In terms of the Schwarzschild radius, i? warp is written 
as (Volonteri et al. 2007) 

^S. = 3.6 x 10 3 a 5/8 ( -^-] A" 1/4 
ifechw V !0 8 M Q J 

x (g) a" 1 ". (15) 

Here a is the Shakura-Sunyaev viscosity parameter and A = 
L/Z/Edd is the Eddington ratio. The accretion luminosity, L, and 
Eddington luminosity, Z/Edd, are defined as 



L = eMc, 
with e denoting the accretion efficiency, and 



(16) 



(10) 



AnGMbhC 46 ( Mbh 

L Edd = = 1.4 x 10 ^— ) erg s , (17) 
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where n ~ 0.3 cm 2 g _1 is the electron scattering opacity. Associ- 
ated to the Eddington luminosity is an accretion rate expressed as 
Mvdd = Lvdd I ec 2 . This is the accretion rate for which the black 
hole radiates at the Eddington luminosity. It is then convenient to 
express the physical accretion rate in units of the Eddington accre- 
tion rate, m = M /M-Bdd- We note that the accretion efficiency in 
MEdd is spin dependent; however, we chose to keep e = 0.1 when 
calculating A/sdd for simplicity. 

Given that the accretion is characterised by a timescale t vi , 
the mass of the disc inside the radius f? warD is 



M d (R warp) — A'lti/^ (^warp)j 

where the accretion timescale is given by 



(18) 



t„ = 



R* 



= 3 x 10° a 



7/8 



f M bh 
V 10 s M G 

xA- 3 / 4 ( ^ 



11/8 



-7/8 



" 3/2 yr. (19) 



We can write the total angular momentum, Jd, passing through 

-Rwarp aS, 

Jd(-Rwarp) ^ A^d(Rwarp)(GAfbh-Rwarp) 1//2 . (20) 

In terms of the anti-alignment criterion this gives. 



Jd 



2J b h aM b 



R* 



Rsc 



1/2 



1 yr 



Rv 



Rs 



1/2 



(21) 



Evaluation of this quantity determines whether or not the anti- 
alignment condition in a misaligned disc is satisfied. 

2.2.3 The case of self-gravity limited discs 

The BH growth process in AGN environments involves vast 
amounts of accreted gas, often comparable to the initial mass of the 
accreting hole. This amount of mass settles onto the accretion disc 
around the hole, which often extends to several thousands of grav- 
itational radii. It is usually assumed that the available mass fuel, 
Mace, is consumed in a single accretion episode, thus providing 
a supply of constant angular momentum (Volonteri et al. 2007). 
In this case, the amount of mass consumed is enough to spin up 
the hole up to a = 0.998, even if the BH initially was maximally 
spinning in a counter-rotating direction. Inevitably, repeated accre- 
tion episodes during major galaxy mergers act to spin up the BH to 
maximum rotation. 

Recently, however. King et al. (2008) argued that the end 
result of the accretion growth channel might be completely dif- 
ferent if we take into account the fact that an accretion disc be- 
comes self-gravitating at some radius, R ss , where its mass exceeds 
M S g ~ (H/R)Mbh- As a result, the mass of the disc is limited by 
its self-gravity to AM cp i so d -C M aC c, which gives rise to a series 
of N ~ Macc/AM cp i so d c well separated accretion episodes. King 
et al. suggest that these accretion episodes are randomly oriented 
around the BH, an assumption supported by observations indicat- 
ing that there is no apparent relation between the accretion disc 
(or radio jet) orientation and the host galaxy disc (Nagar & Wilson 
1999; Kinney et al. 2000). They further argue that this could be the 
result of intense star formation outside R sg randomising the input 
gas direction, which could provide a qualitative explanation for the 
ring of stars seen in the near vicinity of the central BH in the Milky 
Way (Genzel et al. 2003). 



The much smaller angular momentum associated with each 
accretion episode means that, in general, Jd < 2Jbh, so in the 
chaotic accretion model counteralignment occurs in a fraction 
(King et al. 2005) 



/ 



Jd 

2 Jbh 



1/2 



(22) 



of the accretion episodes. Therefore, counter- and co-alignment are 
equally likely outcomes of the Lense - Thirring effect. However, 
accretion of gas in a counter-rotating disc is more efficient in spin- 
ning up the BH since the gas is being dumped onto the BH from 
a larger distance and, thus, carries more angular momentum into 
the BH. Hence, a succession of counter-aligned and co-aligned ac- 
cretion episodes with equal frequency should systematically spin 
down the BH, resulting in a global spin distribution oscillating 
around zero. 

The physics of self-gravitating accretion discs is described by 
Pringle (1981). Briefly, the criterion that the self-gravity of a disc be 
negligible is the requirement that the gravitational force along the 
z direction be dominated by the central BH. This can be expressed 
as the surface density of the disc, E, being negligible compared to 
the quantity MbhH/R 3 , or in terms of the disc mass, 



M d (< R) « ^M bh . 



(23) 



The self-gravity becomes marginally important when Md - 
(H/R)Mbh, which defines the self-gravity radius, _R sg , given by 



Rs 



^ = 1.5xl0 3 e WJ^ 



V 10 s M G 



-26/27 

.,-8/27 14/27 



(24) 



The semi-thickness of the disc in the Shakura-Sunyaev model 
is given by, 



H i Q fi „ 1n -3 -1/5 ( -^bh 



-1/10 



A 



1/5 



R 



1/20 



q- 1/10 . (25) 



-RSchw 

Hence, by replacing R with R sg in expression l !25t . we obtain an 
analytic expression for the disc mass within i? sg , 

H 



AL 



-5- M bh 

H \R=R B g 

( M b 



2.13 x 10 5 e~ 5/27 

23/27 

,5/27 -2/17 



(26) 



V10 8 Mq _ 

As we described in Section l2.2.2l once the disc forms with a non- 
zero misalignment angle about the spin axis of the hole, it will be 
subject to the Lense - Thirring precession. If the radius, R sg , is 
greater than i? W ar P , then 



Jd — Jd ( Rwarp ) 



(27) 



If, however, R sg < i? W ar P , then the entire disc will be subject to 
Lense - Thirring precession, thus aligning itself in the equatorial 
plane of the BH. 



2.2.4 BH-BH binary coalescence 

We now examine the evolution of spin during BH binary coales- 
cence. During the merger of the two binary members, the smaller 
BH adds its spin and orbital angular momentum at the LSO to the 
spin of the larger one. Even if the progenitor holes do not possess 
any orbital angular momentum, the final remnant will preserve the 
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residual orbital momentum of the binary (i.e. the angular momen- 
tum that has not been radiated away). Thus, the remnant will always 
be a Kerr BH. 

Recent breakthroughs in numerical relativity have provided 
robust simulations of BH mergers by solving directly the Einstein 
equations in a fully relativistic frame (Campanelli et al. 2007; Her- 
rmann et al. 2007; Marronetti et al. 2007). These simulations have 
explored the parameter space for different configurations of initial 
masses and spins, allowing accurate measurements of the final spin. 
For example, a merger of two equal mass BHs both with a — in 
a circular orbit results in a Kerr BH with cr ~ 0.69 (Baker et 
al. 2006b). This a robust prediction and is also valid for eccentric 
orbits with eccentricities smaller than 0.4 (Hinder et al. 2008). 

Analytic fits extend the predictions for ar to the entire space of 
parameters and reproduce closely all the available numerical data. 
In the case where the masses are unequal but the spins are zero, the 
final spin can be estimated by the analytic expression 



2V3-, 



- 2.029- 



(28) 



'(1 + q) 2 """" + 

where q — M2/M1 ^ 1 is the mass ratio of the progenitor BHs 
(Berti et al. 2007). 

A more general analytic fitting formula for predicting the final 
spin of any given binary configuration has been provided by Rez- 
zolla et al. (2008a,b). In their analysis, they assume that the final 
spin parameter vector can be expressed as 

f 1 



■(ai + a 2 g + - 



(29) 



where ai, 2 = eJ^/fGMj^) and £ = £'/(MiM 2 ), with £' 
defining the difference between the orbital angular momentum, I, 
when the binary is widely separated, and the angular momentum ra- 
diated away in gravitational waves before the merger, j rad , namely, 

e' = i-3^- OO) 

The vector £' is taken to be parallel to the orbital momentum vec- 
tor throughout the evolution of the binary, an assumption which is 
not strictly valid since the system could radiate away angular mo- 
mentum in a non-symmetric way. However, the error introduced 
by this assumption is relatively small for the binary configurations 
studied in the simulations of Rezzolla et al. A further assumption 
is that the mass radiated away in gravitational waves is negligible, 
as it accounts for only a small fraction (5 — 7%) of the total mass- 
energy of the binary configurations analysed. Under these assump- 
tions, Rezzolla et al. (2008a,b) proposed an analytic expression for 
the magnitude of the final spin, given by 

I = Ti 7 \ 1 l"l a i| 2 + l a 2| 2 <7 4 + 2|a 2 ||ai|g 2 costp 
(1 + q) 1 L 

+2(|ai|cosj9 + \a 2 \q 2 cos £q + l 2 q 2 Y * , (31) 
with the cosine angles ip, 1? and £ defined as 

cos (p = &i • &2, cos 1? = &i ■ £, cos £ = &2 ■ £ . (32) 
The norm of £ is given by, 

W = (J + q 2 ) 2 ^' ai ' 2 + ' a2 ' 2<?4 + 2 l a ill a 2|9 2cos< < 5 ) 
. { Ss/J. + to + 2\ ,. 2 s 

+ 1 i + q 2 J (jail costf + |a 2 |g cos£) 

+2V3 + t 2 ^ + t 3 ^ 2 . (33) 
Here /i expresses the symmetric mass ratio, /i = q/(l + g) 2 = 



M 1 M 2 /{Mi_ + M 2 ) 2 , and the coefficients take the values s 4 = 
-0.129, s 5 = -0.384, t = -2.686, t 2 = -3.454, t 3 = 2.353. 
The final spin as given by Eq. d3U is in good agreement with nu- 
merical data, with residuals of less than 3%. 

2.3 Simulating the spin evolution 

We finally describe our method for calculating the evolution of BH 
spin. During the formation history of the BHs we take into account 
spin changes due to the accretion of gas and mergers, as described 
in Section [2~2l The initial population of seeds is assumed to be 
non-rotating. These seeds grow through accretion of gas during disc 
instabilities, galaxy mergers, and radio mode accretion and through 
BH-BH mergers. 

2.3.1 Thin discs 

To model the physics of the accreted gas, we assume that when 
m ^ 0.01 the gas forms an accretion disc around the BH whose 
physics is described by the standard Shakura-Sunyaev disc model. 
We assume that the gas available to be fed into the BH after a 
galaxy merger or the collapse of a dynamically unstable disc is ac- 
creted over a timescale proportional to the dynamical timescale of 
the galactic bulge that hosts the BH. A fixed proportionality factor 
of 4 is used in this analysis resulting in typical accretion timescales 
of the order of 10 7 — 10 8 yr. Note that, this factor is ~ 10 longer 
than the one chosen by Malbon et al. (2007). The new value is in- 
troduced in order to match the local quasar luminosity function and 
does not affect the galaxy-formation model. The timescale for ac- 
cretion during the radio mode is computed directly by the galaxy 
formation model and depends on the cooling timescale of the gas 
in the hot halo. 

The accretion disc is assumed to be randomly oriented rela- 
tive to the spin axis of the hole. If 8 7^ 0°, 180°, we assume that 
the disc precesses around the BH as described in Section 12.2.21 
In brief, for a given BH mass, Mbh, we determine the numerical 
values of the warp parameters, -R wa r P , Afd(-Rwarp), and the an- 
gular momenta Jd(-Rwarp) and Jbh. The ratio Jd/2Jbh allows us, 
through the criterion in Eq. dl4t . to check if anti-alignment or align- 
ment occurs. We then allow the BH to accrete an amount of gas 
equal to Md(-Rwarp) and evolve the spin using Eq. (|6j. Finally, we 
update the mass of the BH and estimate the hole's new spin and the 
disc precession by calculating the new angles 6 and 9t. This pro- 
cess is repeated until: (a) the disc aligns or counteraligns (6 — 0° 
or 180°) with the hole's spin, in which case we just consume the 
rest of the available gas and evolve the spin according to Eq. ([6j, 
or (b) the accretion disc is entirely consumed, without being able 
to align or counteralign itself with the spin. In almost all cases we 
find that the accretion disc ultimately aligns (or anti-aligns) itself 
on the equatorial plane of the BH. 

In addition to the case where the entire gas reservoir is con- 
sumed within one accretion episode with constant angular momen- 
tum orientation (prolonged accretion), we test the model proposed 
by King et al. (2008) in which the mass of the accretion disc is 
limited by its self-gravity (chaotic accretion). We repeat the same 
steps as before. However, in this case Md ~ M sg . Once M ss is 
consumed, we update the mass of the hole and determine the mass 
of the new disc, ML, computed from the updated BH mass. At this 
point we make the assumption that the disc retains no memory of 
the angular momentum of the initial flow and that it settles into a 
random orientation around the BH with Aid — ML. We then add 
the next M ss until the disc is consumed. 
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It is important to clarify at this point that the two accretion 
models considered here serve as a means for studying the effect 
of accretion on BH spin. The forthcoming analysis that follows 
and results do not rely on any specific model for the physical pro- 
cesses that drive the gas flows towards the vicinity of the BH. Sev- 
eral mechanisms for feeding the SMBH with gas could result to 
similar spin distributions to the chaotic accretion case considered 
here (see for example, Heller et al. 2001 for the properties of non 
self-gravitating gaseous bars in barred disc galaxies and Hobbs et 
al. 2010 for an investigation of the effect of supernova-driven turbu- 
lence on the fuelling of SMBHs in galactic bulges). Here, our main 
goal is to explore the effect of different SMBH spin distributions 
on properties of the host galaxy. 

2.3.2 ADAFs 

Many of the accretion events in our simulations are characterised 
by low, sub-Eddington, accretion rates, with m ^ 0.01. For such 
low accretion rates, the gas flow has low density and is unable to 
cool efficiently since radiative cooling does not balance the energy 
generated by viscosity. Thus, the viscous energy is trapped in the 
gas as entropy and ultimately advected into the hole. This type of 
accretion is known as an ADAF (Rees et al. 1982; Narayan & Yi 
1994; Abramowicz et al. 1995). 

ADAFs have a number of distinct properties that some of them 
will be essential for the analysis in later sections (see Sections 5 and 
6). For example, for an ADAF around a BH, only a fraction of the 
standard accretion luminosity, L — eMc 2 , is emitted as radiation. 
The remainder of the viscously dissipated energy is stored in the 
gas as entropy, resulting in hot flows with almost virial tempera- 
tures. We note that, as shown by Ichimaru (1977), the ions and the 
electrons in an ADAF are not thermally coupled and, thus, reach 
different temperatures. This two-temperature virialised plasma flow 
is optically thin and, for high viscosity parameters (a ~ 0.2 — 0.3), 
it acquires a quasi-spherical geometry around the BH (H ~ R), 
which resembles spherical Bondi accretion. However, despite the 
geometrical similarity, the dynamics of the ADAF are fundamen- 
tally different, since the accretion is entirely due to dissipation via 
viscous forces rather than gravity. 

In general, the accretion flow can have a rather complicated 
structure. For example, it is possible that a thin disc dominates the 
geometry of the outer parts, then switches to an ADAF as the flow 
approaches the BH (Esin et al. 1997). The transition to the thick 
flow depends on the accretion rate and usually occurs closer to the 
BH for higher accretion rates. Such a complex disc model for de- 
scribing the accretion flows during the radio mode is, however, be- 
yond the scope of this paper. Here, we assume a simple configura- 
tion of a quasi-spherical corona around the hole with any transition 
to a thin disc occurring far from the BH. Precession effects due to 
misalignment between the disc and the BH are modelled as in the 
thin disc case. 



2.3.3 Binary mergers 

Finally, we consider the spin changes due to binary coalescence. 
Following the merger of two galaxies in GALFORM, each harbour- 
ing a SMBH, the one hosted by the satellite sinks towards the 
SMBH of the central galaxy and eventually form a binary. In our 
model, BH mergers tend to occur in gas-poor environments, and 
thus, we assume that torques from accreting flows that might be 
present during the formation of the binary do not influence the ori- 
entation of the BH spins prior the merger. We therefore take the BH 
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Figure 3. The distribution of accretion rates for selected BH mass ranges 
at various redshifts. Line styles represent two different accretion modes: 
solid lines for the quasar mode and dashed lines for the radio mode. The 
shaded area represents the regime where the accretion flow is described by 
an ADAF. The accretion rates in the radio mode are truncated at m = 0.01, 
in accordance with our model of AGN feedback. In the present Universe, 
the quasar mode accretion peaks at fa ~ 0.5. However, the peak shifts to 
higher values at higher redshifts, indicating that accretion was more efficient 
in the past. 



spin vectors in Eq. <3U to be randomly oriented relative to each 
other and the orbital angular momentum. In addition, we assume 
that all binaries will ultimately merge on a very short timescale, ef- 
fectively right after the host galaxies merger and induced accretion 
have been completed. 

During the merger, the smaller BH plunges into the larger one 
carrying along its angular momentum at the LSO. At this point we 
neglect any mass loss due to gravitational wave emission, as it cor- 
responds to a small fraction of the total mass-energy of the system. 
In addition, we do not take into account the effects of gravitational 
recoil due to asymmetric emission of gravitational waves in binary 
mergers. In our model, mergers with recoil velocities as high as 
^recoil ~ 500 km s _1 (Libeskind et al. 2006) could displace BH 
remnants. However, we assume that dynamical friction relocates 
them in the galactic centre. Therefore, we do not expect gravita- 
tional recoil to have a significant impact on the spin distribution 
of the remnant BHs. The spin of the final remnant is finally esti- 
mated using the fitting formulae in Eqns. d 3 1 1 > & J331 >, as described 
in Section l2"72"l4l 



3 EVOLUTION OF BH MASS AND SPIN 
3.1 Evolution of BH mass 

We now briefly present our main predictions for BH growth, using 
our updates (see Section 2) to the model developed by Bower et 
al. (2006) and Malbon et al. (2007). In Fig.[3]we show the fraction 
of all BHs that accrete at a given rate in the quasar and radio modes, 
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Figure 4. The correlation of accretion rate with BH mass at redshift zero. 
The median (red-yellow points), 20 — 80 percentiles (dark grey) and 5 — 95 
percentiles (light grey) are shown. The shaded background indicates the re- 
gion were the accretion disc is modelled as an ADAF. Only lower mass 
BHs (< 10 s Mq) i.e. those hosted by lower mass galaxies (typically spi- 
rals) have enough gas left to trigger accretion via a thin disc (rh > 0.01). 
More massive BHs have little cold gas and undergo radio mode accretion 
which (by construction) has rh < 0.01. 



Figure 5. Contribution to the final BH mass at different redshifts from 
gas accretion in quasar mode (solid green line), radio mode (dashed blue 
line) and BH-BH mergers (dashed-dotted red line) in the updated Bower et 
al. (2006) model. The lines represent the medians and the shaded region the 
20-80 percentiles of the distributions. The dominant channel of BH growth 
for BHs with Mbh $ 10 8 Mq is the quasar mode. Above that mass, merg- 
ers and radio mode become the dominant growth channels. 



at different redshifts. All BHs that accrete in the radio mode dur- 
ing the redshift bin are included, but the quasar mode is episodic 
so we only include objects which experience a starburst within a 
given redshift bin. Our updated parameters give a longer duration 
for the accretion episode in the quasar mode (4 times rather than 
0.5 times the dynamical timescale of the bulge) than assumed in 
Malbon et al. 2007. This results in the majority of BHs with mass 
10 6 — 10 s Mq which accrete in the quasar mode at z = having 
rh ~ 0.5, as observed (Heckman et al. 2004). At higher redshift 
the distributions shift to higher mass accretion rates, with a mean 
rh ~ 3 at z = 6. This implies that in a hierarchical universe, accre- 
tion of gas was more efficient at early epochs, so the rate of growth 
of BHs is faster in the past. By contrast, accretion activity during 
the radio mode peaks at very low accretion rates (rh ~ 10 -4 ) at 
z = 0, with a long tail to lower values together with some spread 
to the maximum allowed value of 0.01 (imposed by the AGN feed- 
back model). This upper limit to the mass accretion rate limits the 
shift to higher rh at higher z, so the main evolutionary trend is that 
there are fewer very low rh radio mode objects at high z. 

This model results in a bimodal mass accretion rate distribu- 
tion for all BH masses and redshifts. Activity triggered by accretion 
of cold gas from the galaxy disc, replenished by galaxy mergers, 
has a mean mass accretion rate which is ~ 1000 times higher than 
that of activity fed by radio mode accretion. However, the relative 
importance of these two modes varies with BH mass. This is shown 
in more detail in Fig. [4] which shows the 20 — 80 and 5 — 95 per- 
centile mass accretion rates as a function of BH mass at z = 0. 
More than 5 per cent of the lower mass AGN have accretion via 
a thin disc (rh > 0.01). However, this fraction drops with mass 
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Figure 6. The Mbh — Mb u i gc relation predicted by the updated Bower et 
al. (2006) model (solid line). The shaded areas indicate the 10 — 90 (light) 
and 20 — 80 (dark) percentile spread of the theoretical predictions. The 
observational data is taken from Haring & Rix (2004). 
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Figure 7. The predicted mass functions at z = 0, 2, 4 and 6 for the BH 
population in our simulation. 

at BH masses of ~ 10 8 M Q . More massive galaxies are hosted 
by large ellipticals which are gas poor and are dominated by radio 
mode accretion. 

Fig- shows the contribution from the different growth chan- 
nels to the final BH mass at different redshifts. At high redshifts, 
SMBHs with masses up to ~ 10 8 Mq build their mass almost 
exclusively through the accretion of cold gas (quasar mode). This 
forms either a thin disc or an ADAF depending on m (see also 
Fig-EJ. Accretion of gas from the hot halo of massive galaxies dur- 
ing the radio mode always occurs with m ^ 0.01 by construction, 
thus giving rise to an ADAF. This low upper limit to the mass accre- 
tion rate in this mode means that it makes comparatively little im- 
pact on BH growth, except for masses greater than 10 8 Mq. Note 
that in the original Bower et al. (2006) model, the radio mode is re- 
sponsible for building most of the BH mass above 10 s Mq, while 
in the updated version presented here, both BH mergers and radio 
mode contribute significantly for high mass BHs. 

In Fig. [6] we compare the relationship between the BH mass 
and host bulge mass predicted by our model to the observational 
data. The model predicts an almost linear relation and reproduces 
the observations reasonably well. The scatter in BH mass is ~ 
1 dex for bulge masses below 10 11 Mq and becomes gradually 
smaller for higher bulge masses. The resulting BH mass functions 
are shown in Fig. [7] These show that the present day Universe is 
dominated by BHs with masses of ~ 10 7 - 10 s M which is in 
good agreement with the recent results from the Sloan Digital Sky 
Survey (SDSS) of Heckman et al. (2004). The shape of the mass 
functions are fairly similar at all redshifts. The main difference oc- 
curs at the high-mass end where the BH mass function extends to 
higher masses at low redshifts. For Mbh ^ 10 Mq the mass func- 
tion is almost constant and flat up to z — 2. This indicates that 
these BH are already in place at z = 2. Above that mass, there is a 
sharp decrease in the space density of BH at z = 2 and z = 4. At 
z = 0, the mass function has a shallower slope. The estimated BH 
mass density at z = is pbh = 3.84 x 10 s M Q Mpc~ 3 , which is 



in good agreement with the value implied by the X-ray background 
(Fabian & Iwasawa 1999). 

3.2 Global BH spin distributions 

The histograms in Fig. [8] show the distributions of BH spins at 
different redshifts (left panel) and for different mass ranges (right 
panel) for all galaxies that host BHs with Mbh > 10 6 M Q . The 
top and bottom panels correspond to the chaotic and prolonged ac- 
cretion models respectively. The BH spin distributions predicted 
by the two accretion models are already established at z = 6 and 
do not change appreciably with time. The major channel for BH 
growth at z = 6 is the quasar mode and this generally takes place 
via a thin disc. The accreted mass in each episode is generally larger 
than the mass of the BH, so this spins it up to maximal in the pro- 
longed accretion model. The break up of the accretion into much 
smaller individual events in the chaotic model leads to a random 
walk spin distribution around a low spin value. Nonetheless, the 
plots in Fig. [5] show that BH-BH mergers contribute to the most 
massive BHs at late redshifts, as does accretion with m < 0.01 
in both the radio and quasar mode. In order to gain further insight 
into the effect of these growth channels we show the different spin 
distributions as a function of BH mass at z = in the right hand 
panels. 

As already mentioned, low mass BHs are built via thin-disc 
accretion in the quasar mode (Fig.0 with mean rh ~ 0.5. For a 
BH of mass 10 6 Mq a typical accretion event in the quasar mode 
contains 10 6 Mq of gas, so in the prolonged accretion model this 
will spin the BH up to maximal. However, in the chaotic accretion 
model in which the size of the disc is limited by its self-gravity, 
the mass of the gas that settles on the disc cannot be more than 
~ (H/R) x 10 6 M ~ 10* M Q . Thus, a 10 6 M© accretion 
event consist of ~ 100 accretion episodes with AM op j S odc <C 
Mbh, all randomly oriented around the BH. As a consequence, 
10 6 — 10 s Mq BHs experience a net spin-down to modest values 
centred around a ~ 0.15. 

The growth channel for more massive BHs (Mbh > 10 8 Mq) 
includes a large fraction of objects which accrete in the radio mode 
or quasar mode with m < 0.01 i.e. via a geometrically thick 
ADAF. However, we do not include full modelling of this pro- 
cess as our spin evolution calculations always assume a Shakura- 
Sunyaev thin disc rather than the appropriate ADAF equations. We 
do not expect this to have a large effect on the spin distributions as 
most of the BH mass and spin are built up from high mass accretion 
rates. 

The mergers between massive BHs have a significant impact 
on the final spin distributions. According to the analytic fitting for- 
mula, Eqn. (EB from Rezzolla (2008, 2009), the post -merger spin 
of the final remnant depends strongly on the masses of the pro- 
genitors. For example, BHs acquiring a final spin greater than 0.69 
after a merger are the end product of binaries of comparable mass 
(q ~ 1) in which the members already had significant spins (re- 
call that a f = 0.69 is the final spin of a binary of equal-mass 
non-spinning BHs). Mergers between BHs of comparable mass are 
common only for the most massive BHs (Mbh > 10 Mq) in our 
simulation^ This is because these BHs are hosted by massive el- 

3 At redshift zero nearly 19% of the BHs with M bh 10 6 M Q have 
experienced at least one merger. The low merger rate of BHs merely re- 
flects the relative minor effect that mergers have in the formation of galac- 
tic spheroids (except in very massive ellipticals) in the Bower et al. (2006) 
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Figure 8. Model predictions for the final BH spin distributions, (a) The spin evolution due to BH mergers and accretion and (b) the distribution of BH spins 
at z = in different mass bins in the chaotic (upper panels) and prolonged (lower panels) accretion models. The main difference between the two accretion 
modes is that the chaotic case results in low mass BHs having spins distributed around zero, whereas the prolonged case results in a distribution peaked at 
around a f = 0.998. 



liptical galaxies that have experienced at least one major merger 
(merger between galaxies of equal mass) in their past history (Parry 
et al. 2009). Hence, their BHs, which are of similar mass, will ac- 
quire a final spin close 0.69 or higher (depending on the value of 
the initial spins) after they coalesce. In contrast to the most mas- 
sive BHs, lower mass BHs (Mbh < 10 s M ) which experience a 
merger, are involved mostly in minor mergers (q <C 1). These do 



model, as found by Pany et al. (2009). For example, at redshift zero about 
6% of the BHs with Af^h 5* 10 6 Mq have experienced more than one 
merger and only ~ 8% have had a recent major merger 



not have a significant impact on the final spin of the remnant which 
is dominated by the more massive member of the binary. 

When the effects on BH spin of mergers and accretion are 
combined, we find the following result. BHs of mass > 5 x 10 s Mq 
that have experienced a major merger will acquire spins greater 
than ~ 0.69; since these mergers often occur in gas poor environ- 
ments, the spin distribution for these masses is not significantly al- 
tered by the accretion. For lower mass BHs, accretion quickly erases 
the characteristic post-merger spins, because the growth is domi- 
nated by intense accretion during the quasar mode. For chaotic ac- 
cretion, this results in a bimodal distribution of spins for BHs with 



© 0000 RAS, MNRAS 000, 000-000 



12 Fanidakis et al. 




Figure 9. The correlation of spin with BH mass (a) for the prolonged and (b) the chaotic accretion models. The solid lines show the medians and the shaded 
areas the 10 — 90 (light) and 20 — 80 (dark) percentile spreads of the distributions. 



mass 10 8 — 10 10 Mq . The first peak is located at a ~ 0.15 and cor- 
responds to BHs that have had a major merger accompanied by an 
accretion episode in which the post merger spin is quickly erased 
and kept low. The second peak is located at a f ~ 0.7 — 0.8 and 
corresponds to BHs that experienced a merger in the absence of ac- 
creting flows and acquired a final spin characteristic of post-merger 
BH remnants. 

We note that BHs in our cosmological model have low merger 
rates at high redshift (z ~ 6). As redshift decreases, the merger 
rate increases and reaches a maximum at z ~ 1.5. Eventually, at 
2 = nearly 19% of the BHs with M hh ^ 10 6 M Q have expe- 
rienced at least one merger. Therefore, the significance of mergers 
increases with decreasing redshift. This accounts for the decrease 
with redshift of the rapidly rotating BHs in the prolonged model. 

We finally show in Fig.|9]the correlation between the BH mass 
and spin. This figure confirms the conclusion of Fig. [8] BHs that 
grow their mass through chaotic accretion display a strong correla- 
tion between their mass and spin. In this case we find that the BH 
mass correlates with the host-galaxy morphology such that small 
BHs are usually found in spiral galaxies, whereas massive BHs are 
found in massive ellipticals. Thus, the apparent correlation between 
BH mass and spin indicates a strong correlation between BH spin 
and host-galaxy morphology. As a consequence, we expect to find 
rapidly rotating BHs in principle at the centres of very massive el- 
liptical galaxies. However, these are also the objects which have 
low mass accretion rates (see Fig.|4j», and so predominantly accrete 
via an ADAF whereas some fraction of the lower mass, low spin 
BH can accrete via a thin disc. By contrast, in the case of prolonged 
accretion there is no apparent correlation between the BH mass and 
spin since most of the BHs exhibit very high spin independently of 
mass. Thus, rapidly rotating BHs are found in most types of galaxy. 
The correlation of mass and mass accretion rate (Fig.[4} remains the 
same, so the most massive ellipticals harbour the most massive BHs 



and accrete via an ADAF while a small fraction of the lower mass 
BHs in spirals can accrete via a thin disc. Note that the correlation 
between Mhh and spin for BH masses greater than 10 9 Mq is sim- 
ilar to that predicted by the chaotic model since these BHs include 
more growth through BH merger activity, which is independent of 
the accretion model. 



4 OPTICAL LUMINOSITIES OF AGN AND THE 
QUASAR LUMINOSITY FUNCTION 

We now explore whether our models are capable of explaining the 
observed properties of AGN. In particular, we present our predic- 
tions for the quasar luminosity function assuming that quasar ac- 
tivity is driven by accretion onto a SMBH. In our model, SMBHs 
with m > 0.01 are assumed to accrete through a thin disc, whereas, 
those with accretion rates below this are assumed to accrete through 
an ADAF. In the thin-disc regime the bolometric luminosity is 
given by the standard expression, 

L = eMc. (34) 

The accretion efficiency, e, is taken to vary with the hole's spin 
according to Eq. y). When the accretion rate becomes super- 
Eddington (m > 1), the bolometric luminosity is limited to 
(1 + lnm)Z/Edd (Shakura & Sunyaev 1973). However, we do not 
restrict the accretion rate if the flow becomes super-Eddington. For 
the ADAFs, we assume that the plasma acquires a two-temperature 
configuration. The bolometric luminosity of the disc in this case is 
given by Mahadevan (1997), 

'--"-"(I) (S) {&) *> 
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Figure 10. The model quasar luminosity function at z = 0.5 for the pro- 
longed (dashed-dotted line) and chaotic (solid line) accretion models. The 
luminosity function includes only AGN assumed to be powered by a thin 
disc (m > 0.01) that accrete in the quasar mode. The best fits to the obser- 
vational data require that 8% (prolonged model) and 20% (chaotic model) 
of the bolometric luminosity be emitted in the bj-band. It further requires 
that 40% of the total number of AGN be obscured in both accretion models. 
The quasar luminosity functions from the 2dF (Groom et al. 2004) and 2dF- 
SDSS QSO survey (Croom et al. 2009) over the range 0.4 < z < 0.68 are 
shown with circles and boxes respectively. 



where, /3 is related to the Shakura-Sunyaev viscosity parameter a 
through the relation a w 0.55(1 - /3) (Hawley et al. 1995). The 
value of a is taken to be to 0. 1 for all objects in our samples. 

Observationally, the quasar bolometric luminosity function 
can be estimated by combining measurements in many different 
wavelengths (Hopkins et al. 2007; Shankar et al. 2009). We can 
compare the predictions of our model to the estimated bolometric 
luminosity function simply by counting the number of quasars ra- 
diating at a given luminosity bin. In this analysis, however, we will 
attempt to calculate the optical 6,j-band luminosity function since 
we are ultimately interested in comparing the strength of the radio 
jet emission of an AGN relative to its disk optical emission (see 
Section 6). 

The conversion factor between bolometric accretion luminos- 
ity and &j-band luminosity, f bj , is assumed to be constant and its 
value is adjusted to match the observed foj-band luminosity func- 
tion. We take f b , = 0.08 for the prolonged accretion and f bj =0.2 
for the chaotic accretion model. The substantial difference between 
the conversion factors for the two accretion models reflects the 
strong dependence of the bolometric luminosity on the spin through 
the accretion efficiency. An increase in spin means that the disc 
can extend closer in towards the BH, so the same mass accretion 
rate produces more UV emission for the same 6,j-band luminosity. 
Hence the b j -band carries a smaller fraction of the bolometric lu- 
minosity. We check our conversion factors derived from fitting the 
quasar luminosity function to those derived from a Novikov-Thorne 
(fully relativistic) version of the Shakura-Sunyaev disc equations, 
assuming that the bj band is centred at 4400 A, with a FWHM of 



980 A. We find that our assumption is in good agreement with the 
theoretical calculations as long as we assume that the disc luminos- 
ity is Eddington limited and similarly they are in good agreement 
with the values of 0.07—0.15 suggested by Elvis et al. (1994) from 
analysis of the SEDs of a large sample of quasars. Thus, we adopt 
these values for the rest of our calculations. 

In deriving the luminosity function, we assume that only thin- 
disc accreting objects have the intense ionising flux which leads to 
their identification as quasars. We further assume that a fraction of 
40% of the total number of quasars in obscured in the &j-band by 
the torus (Polletta et al. 2008). Finally, absolute 6,j -magnitudes in 
the Vega system are obtained using the relation, 

M bj = -10.44 - 2.51og(L 6j /10 40 ergs~ 1 ). (36) 

Fig. [TU] shows the resulting quasar luminosity function at red- 
shift zero for both the prolonged-accretion model, and the chaotic- 
accretion model. The models are compared to the observed lumi- 
nosity functions in the redshift range 0.4 < z < 0.68 estimated 
from the 2dF (Croom et al. 2004) and 2dF-SDSS QSO survey 
(Croom et al. 2009). The two models agree with each other, al- 
though the prolonged-accretion model predicts a somewhat higher 
volume density of quasars in the —21 ^ Afb T —23 and 
-24 < M bj < -25 range. The inflection at M bl ~ -24 that 
separates the two regions is a reflection of the 10 s M© dip in the 
Mbh — a correlation seen in Fig. [9] Both models predict steeper 
faint end slope (—20 ^ M bl —22) than the observed luminosity 
functions. A similarly steep faint end is also seen in the soft X-ray 
luminosity function at z ~ 0.4 (Miyaji et al. 2001; Hasinger et 
al. 2005). The differences between the two models in this regime 
arises because the BHs in the prolonged model that contribute to 
the faint end have much higher spins, and thus higher bolometric 
luminosities than in the chaotic model. 

In conclusion, we can reproduce the optical luminosity func- 
tion reasonably well with either the chaotic or prolonged accretion 
models. A consistency check is provided by calculating the aver- 
age accretion efficiency. Assuming only the accreting sources con- 
tributing to the quasar luminosity function we find (e) = 14.2% 
((a) = 0.86) for the prolonged model and (e) = 6.1% ((a) = 
0.12) for the chaotic model. The latter is in good agreement with 
(e) = 6.7% inferred by optically selected AGN (Martfnez-Sansigre 
& Taylor 2009) and consistent with the low efficiency that charac- 
terises accretion at z — suggested by Wang et al. (2009). 



5 AGN RADIO LOUDNESS AND THE SPIN PARADIGM 

In this Section we investigate how the total radio output of an AGN 
is related to the accretion process and the central BH properties, 
and we present a model where we combine standard accretion disc 
theory (including ADAFs), BH evolution and radio jet production 
to understand the observed AGN radio loudness. 

5.1 Modelling the jet emission in AGN 

We construct a model for studying the radio loudness of AGN based 
on our current understanding of the formation and acceleration of 
jets. The starting point is the galaxy formation model described 
in Section|2] Galaxies become active every time the central SMBH 
experiences an accretion episode, allowing us to track the evolution 
of BH mass, spin and mass accretion rate distributions as described 
in the previous section. We combine this with the prescriptions for 
jet luminosity in the BZ model (Meier 2002). For consistency, we 
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also investigate the more general case where the jet taps energy 
from both the BH and the accretion disc (Blanford & Payne 1982), 
which is described in the hybrid model presented by Meier (2001). 

5.1.1 The BZ jet model 

AGN jets in the BZ jet model are exclusively powered by extraction 
of the rotational energy of the BH and, thus, they are formed only 
in AGN hosting rotating SMBHs. The mechanical energy of a jet 
is proportional to the square of the poloidal magnetic field, B po \, at 
the horizon of the BH (Blandford & Znajek 1977), 

L\ et OC Bp 0l A/ b > 2 . (37) 

This is a second-order perturbative solution for the spin parame- 
ter a and is used to approximate the solution for slowly rotating 
BHs, a <C 0.998. By including higher order corrections it can be 
shown that the dependance of Lj e t on a can be much steeper for 
rapidly rotating BHs (Ljct ~ a 4 ; Tanabe & Nagataki 2008; see 
also Tchekhovskoy et al. 2010 for a recent study of the dynamical 
range of Lj ct predicted by these solutions). For this analysis we will 
adopt the original approximation of Blandford & Znajek (1977) for 
simplicity. 

There have been several models that give the strength of B po \ 
around a rotating BH. In most cases, B po i is expressed in terms of 
the azimuthal component, B^, as B po \ ~ (H / ' R)B^. Therefore, 
the strength of the jet depends critically on whether the disc is geo- 
metrically thin or thick. For ADAFs {H ~ R) and thin discs (TD; 
H -C R), this yields (see Meier 2002 and references therein): 

i,..^ = 2 x 10- (^)(^)a» «,.-, (38, 

In the super-Eddington regime (m ^ 1), we assume that the 
flow remains in a thin disc state as there are as yet no models to de- 
scribe the behaviour of the radio jet in this regime. The value of the 
viscosity parameter, a, in the ADAF and thin-disc regimes is set to 
0.1. Note that, according to Eqn. l |38t the mechanical jet power at 
the top of the ADAF branch is Lj ct = a 2 Lboi ~ 0.01a 2 Le<m ergs 
s _1 . This upper limit to the jet luminosity in the ADAF regime is 
the origin of our revision to the fraction of the Eddington luminos- 
ity available for jet feedback into the halo (see Section 2). For low 
spin BHs, the BZ jet power drops substantially below this. 

5.1.2 The hybrid jet model 

In the BZ mechanism, the jet is assumed to be powered directly by 
the extraction of the rotational energy of the BH, neglecting any 
energy that could be extracted from the disc. However, the field 
lines frozen to the accreting matter on the disc may generate col- 
limated outflows as a response to the differential rotation of the 
plasma, even in the case when the BH is not rotating (Blandford & 
Payne 1982). The outflow in this case is powered by the extraction 
of the rotational energy of the disc, rather than the BH. Thus, in 
the general case, both the disc and the hole could contribute to the 
production of jets. 

Based on these considerations, Meier (2001) proposed a hy- 
brid model for the radio loudness of AGN in which the jet lumi- 
nosity depends on the rotation of both the BH and the accretion 
disc. The model includes two distinct accretion states, the ADAF 



and thin disc, in which the BH can be rotating or non-rotating. The 
hybrid model has a weaker dependence on the spin of the BH than 
the standard BZ mechanism. When a / the jet luminosity scales 
with spin as 

£5,A r DAF oc (0.55/ 2 + 1.5/a + a 2 ), (40) 

£j2,TD oc (1 + 1.1a + 0.29a 2 ). (41) 

where / and g are parameter related to the angular velocity of the 
disc and the azimuthal magnetic field respectively (see Meier 2001 
for more details; also Nemmen et al. 2007 for a fully relativistic 
account of the hybrid model). The model appears to be a more re- 
alistic interpretation of the jets seen from MHD simulations, since 
these are still able to power jets even from Schwarzschild BH accre- 
tion. We have tested this model by coupling it to our predicted spin 
distributions. However, we find that the weak dependence on spin in 
this model results in small dynamic range in radio luminosity, and 
the results are very similar to the predictions of the BZ jet model for 
the prolonged accretion case (see Section|6j. We therefore prefer to 
contrast the prolonged accretion model with the chaotic accretion 
model coupled to the spin dependence of the BZ jet model. 

5.2 From jet power to radio luminosity 

The mechanical BZ jet luminosities presented in the previous sec- 
tion are converted into radio luminosities according to the theo- 
retical calculations of Heinz & Sunyaev (2003) (see also Falcke 
et al. 1995). Heinz & Sunyaev derived the non-linear dependence 
between the jet flux and the physical parameters Afbh and rn for 
radiatively efficient and inefficient flows. These authors found that 
the core flux at frequency, F v , of a jet scales as M^irh^ 2 , with 
£i = £ 2 = 17/12 for ADAF systems and £i = 17/12, £2 = 
when the disc is radiation-pressure supported A similar depen- 
dence on the accretion rate and BH mass seems to be implied by 
the fundamental plane of BH activity in the case of radiatively inef- 
ficient discs (Merloni et al. 2003; Falcke et al. 2004). Hence, when 
jets are launched in ADAFs we assume that 

£r,adaf oc (M hh m) 1A2 . (42) 

When the ADAF collapses to a thin disc we assume that for the 
range of accretion rates we are interested in, 0.01 ^ m < 100, the 
flow is radiation-pressure supported and thus, 

L r ,td oc A/^ 42 (43) 

Since we know how Lj e t depends on Afbh and rii we can work 
out the relation between Lr and Lj ot for both accretion regimes. 
For example, in the BZ model, Eqns. l |38t and l |39t in combination 
with Eqns. d42t and d43T > give 

£r,adaf = Ai(Afbhm)°' 42 Lj ot ,ADAF, (44) 

£r,td = A 2 M°h 2 m~ 12 Lj ct! TD- (45) 

where Ai and A2 are normalisation constants. In order to con- 
strain the values of these constants, we assume that A1/A2 — 100 
which arises from the fact that the mechanical power of a jet in 
an ADAF is approximately 100 times higher than that of a jet in 
a thin disc (Eqns. [38] & [39} . Ai is then constrained by fitting the 
predictions of each accretion model for the bright end of the radio 

4 The values for the scaling indices §1 and §2 are calculated for a flat spec- 
trum. 
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Figure 11. The predictions of the BZ model for the optical and radio lumi- 
nosities of an AGN. Results are shown for an AGN hosting a SMBH with 
M hh = 10 8 M Q and rotating with spin a = 0.01,0.1,0.5 and 0.998. 
The viscosity parameter a is fixed to 0. 1 in both the ADAF and thin-disc 
regimes. We represent the different accretion regimes with different back- 
ground shadings. The sharp drop in radio luminosity seen in all cases is due 
to the transition form the ADAF to the thin disc regimes when the accre- 
tion rate exceeds O.OlLrjdd. The rightmost shaded area depicts the region 
where accretion becomes super-Eddington. In this region, the optical lumi- 
nosity is set to (1 + In m)L Edd . The discontinuity above 10 42 erg s _1 in 
accretion luminosity is due to the fact that thin discs are radiatfvely efficient 
and thus, for the same accretion rate, are more luminous than an ADAF. 

luminosity function to the observations (see Section 15.41 ). When 
we consider the prolonged accretion model this gives A\ = 0.05 
and A2 = ^i/IOO = 0.0005. For the chaotic accretion model, 
Ai = 0.07 and A 2 = At/100 = 0.0007. The uncertainty in the 
parameters Ai and A2 introduces an arbitrariness into the model, 
which, however, is unavoidable since the normalisation of L R is 
not calculable from first principles. 

5.3 The log L B - log L R plane for a 10 8 M Q BH 

The theoretical predictions of the BZ model for the evolution of 
AGN jets on the log Lb — log Lr, plane are illustrated in Fig. II II 
The optical luminosity, Lb, is obtained as discussed in Section 4. 
The jet strength is explicitly determined by the three fundamental 
parameters of our model, the spin, the accretion rate and the BH 
mass (and also by the values of the normalisation constants Ai and 
A2). The predictions of the models are shown for a BH mass of 
10 s M spinning at a = 0.01, 0.1, 0.5 and 0.998. The different ac- 
cretion regimes are represented by different background shadings. 
As shown in the diagram, a BH that accretes through an ADAF 
develops jets whose radio luminosity increases with the accretion 
rate on the log Lb —log Lr plane. When the accretion rate becomes 
higher than 0.01, the accretion flow enters the thin-disc regime and 
the jet collapses by a factor roughly equal to the difference in H/R 
between the two very different accretion flows. 

In the thin disc regime, the mechanical jet power increases 
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Figure 12. The radio luminosity function of AGN at 1 .4 GHz. The lines 
show our model predictions for the prolonged (dashed-dotted) and chaotic 
(solid) accretion models. The symbols show observational measurements 
for the local radio luminosity function from Sadler et al. (2002) for the 2dF- 
GRS (circles) and Best et al. (2005) for the SDSS (boxes). The predictions 
of the BZ model, in combination with the chaotic accretion model, repro- 
duce the observed radio luminosity function reasonably well. 

with TTi but this does not translate into an increase in observed ra- 
dio flux. Instead, the radio flux remains constant, while the optical 
flux increases with mass accretion rate. Thus, the much lower radio 
emission is sustained throughout the thin disc regime and we do not 
expect objects powered by thin discs to be radio loud. Our model 
extends this behaviour to super-Eddington mass accretion rates, al- 
though we note that this is almost certainly an underestimate of the 
jet and radio luminosities as H/R for such flows increases again 
to ~ 1. Thus, there is a clear accretion mode switch in the models. 
The brightest radio sources are the highest spin BHs accreting at 
the maximum rate for an ADAF. Since this is L/L^dd — 0.01, this 
maximum radio luminosity scales also with the mass of the BH. 



5.4 The radio luminosity function 

A first test of the jet model is the radio luminosity function of AGN. 
In Fig. [12] we show the predictions of the BZ jet model for the ra- 
dio luminosity function in both the prolonged and chaotic accre- 
tion cases and compare with the estimates from the Two-degree- 
Field Galaxy Redshift Survey (2dFGRS; Sadler et al. 2002) and 
from the SDSS (Best et al. 2005]3 The model radio luminos- 
ity function is derived for systems powered by ADAFs and thin 
discs, including super-Eddington objects, in a simulation volume 
of V = 1.25 x 10 8 Mpc~ 3 . 

The coupling between the different spin distributions and the 
BZ mechanism gives interestingly different predictions for the two 



5 The results from 2dFGRS have been adjusted to the cosmology adopted 
by Best et al. (2005). 
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accretion models. The chaotic model shows a fairly good agree- 
ment with the observations throughout the entire energy spectrum. 
The faint end is dominated mainly by galaxies that host slowly ro- 
tating BHs whereas the bright end is dominated by galaxies that 
host rapidly rotating massive BHs. The majority of these BHs ac- 
crete in the radio mode. A characteristic break in the slope at 
~ 10 25 W Hz -1 separates the two regimes. 

In contrast, in the prolonged accretion case the same sample of 
galaxies shows a different distribution of luminosities. The model 
overproduces the number density of faint radio sources by a factor 
of ~ 3 compared to the chaotic model, due to the fact that these lu- 
minosity bins are populated by galaxies with rapidly rotating BHs 
whose radio luminosity is too larger. This trend extends down to lu- 
minosities of 10 25 W Hz -1 . Above that luminosity, the model is in 
reasonably good agreement with the observations. Since the mech- 
anism that spins those BHs up, namely the BH mergers, is present 
in both models, the two predicted luminosity functions agree with 
each other very well at 10 25 - 10 26 W Hz -1 . Note that the param- 
eters Ai and A2 introduced in Section |5?2l affect only the normal- 
isation of the radio luminosity function. The shape is a prediction 
of the model and depends only on the distribution of BH masses, 
spins, and accretion rates. 



6 PREDICTIONS FOR THE L B - L R AGN ACTIVITY 

6.1 The distribution of galaxies on the optical-radio plane 

We now explore the distribution of AGN on the log Lb — log Lr, 
plane. We consider a sample of galaxies brighter than L* (My — 
—20.5). This distribution is displayed in Fig.[T3] where we plot the 
L-band versus the radio luminosity of AGN for both the prolonged 
and chaotic accretion models (left and right plots respectively). All 
galaxies are volume weighted and those accreting through a thin 
disc are weighted according to the fraction of their lifetime dur- 
ing which they are active. The colour coding corresponds to dif- 
ferent types of accretion: blue shows AGN powered by an ADAF, 
red those powered by a thin disc and green indicates the super- 
Eddington objects. 

Along with the model predictions, we plot the radio and opti- 
cal luminosities of a sample of AGN-powered radio sources stud- 
ied by Sikora et al. (2007). We include only local objects i.e. those 
at z < 0.14, so this excludes the very brightest radio and optical 
sources. The evolution of the sources across cosmic time will be 
explored in a future paper. 

The data form two distinct sequences on the log Lb — log Lr 
plane. The upper sequence represents radio-selected galaxies and 
the lower sequence optically-selected objects. The objects in the 
upper sequence are on average ~ 3 orders of magnitudes more ra- 
dio loud than those in the lower sequence. These radio bright ob- 
jects are Fanaroff and Riley class I (FR-I) galaxies, broad line re- 
gion galaxies (BLRGs) and radio loud quasars (RLQs). These are 
thought to be powered by very massive BHs (Mclure & Dunlop 
2002), hosted by elliptical galaxies. Objects in the lower, radio- 
quiet sequence include Seyferts, low ionisation emission regions 
(LINERs) and radio-quiet quasars (RQQs). 

The observational data suffer from multiple selection effects, 
which may be responsible for the apparent dichotomy seen in the 
distribution (see discussion by Sikora et al. 2007). The fact that 
some FRI galaxies have much higher radio power than predicted 
by our model could be due to the fact that these data include lobe 
power, whereas we only model the core radio emission. However, 



the locus of the data points sets upper limits for the radio luminosity 
of radio-quiet and radio-loud objects, and provides a rough visual 
description of the general bulk properties of the different popula- 
tions. 

In both accretion models, AGN span a wide range of optical 
and radio luminosities. AGN powered by ADAFs generically pro- 
duce more radio and less optical emission than those powered by 
a thin disc. Thus, the accretion mode switch at m = 0.01 from 
ADAF to thin disc produces a marked transition in radio loudness. 
All our models show that the radio-loud sequence lies at the upper 
end of the radio luminosity envelope for SMBH accreting via a hot 
flow. Conversely, the radio-quiet locus lies predominantly on the 
thin disc points, although the lowest luminosity radio-quiet AGN 
(LINERS) lie on the lower end of the ADAF regime. Thus, the ma- 
jor switch between radio-loud and radio-quiet is the change in jet 
properties as the accretion flow collapses from a hot, geometrically 
thick configuration, to a cool, geometrically thin disc (Jester 2005). 

Nonetheless, there can also be an additional spin dependence 
which enhances the difference in radio properties, although this de- 
pends strongly on the details of the jet modelling and the spin dis- 
tribution. The BZ jet models are strongly dependent on spin, so 
their radio luminosity amplifies differences in the spin distribution. 
In the prolonged accretion model almost all SMBHs have maximal 
spin, so there is little dispersion in radio flux. Conversely, for the 
chaotic model, the correlation between mass and spin means that 
the low mass and hence lower optical luminosity objects have dra- 
matically lower radio power than the higher mass BH in both the 
ADAF and thin disc regime. 

There are also subtle differences in optical luminosity between 
the prolonged and chaotic models, as the higher efficiency of high 
spin accretion means that the same mass accretion rate gives rise to 
a higher luminosity.Thus, there are more super-Eddington sources 
in the prolonged accretion model, in which low mass objects have 
high spin, than in the chaotic model. 

6.2 The distribution of galaxies on the 7L — A plane 

In addition to the AGN optical and radio output, the model allows 
us to study how the radio loudness 3L of AGN depends on vari- 
ous physical parameters such as the BH mass and the Eddington 
ratio A = Lboi/LEdd- The radio loudness measures the radio to 
optical flux ratio, 3? = L„ R /L„ opt (Richards et al. 2006). Fol- 
lowing the definition of 3? in Sikora et al. (2007), we consider 
the total radio flux at 5GHz, Lr = vs Gllz .L V5 > an d therefore 
we express the radio loudness as 31 = (^b/K5 GHz ) X Lr/Lb = 
1.36 x 10 Lr/Lb, given that the B-band is centred at a wavelength 
of 4400A. 

We plot the theoretical predictions for the radio loudness of the 
AGN in our sample along with the observational data set of Sikora 
et al. (2007) in Fig. [14] Our predictions suggest a clear inverse cor- 
relation between 3i and A with a substantial scatter in both accretion 
models (upper panels in Fig. 1 141). In the ADAF regime, the distri- 
bution of objects shows a correlation of the form 3? oc A -0 ' 4 (one 
can derive the correlation between 31 and A using Eqns. t44t . d 3 81 > 
and l |35t ) which is driven mainly by the strong dependence of the 
disc luminosity on the accretion rate, Lboi oc m 2 . The correlation 
between 3c and A becomes steeper in the thin-disc regime, where 
the data approximately follow 31 oc A -1 . The dependence of 3c on 
A in the thin-disc regime arises from Ji oc L^ ol oc mT x oc A -1 . 

In both models, the radio-quiet AGN (log 3c < 1, following 
Ho 2002) are preferentially found in the log A > — 2 regime, while 
the radio-loud sources (log 3J > 1) populate exclusively the log A < 
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Figure 13. Radio luminosity vs. /3-band nuclear luminosity for the BZ model in (a) the prolonged and (b) the chaotic accretion models. Only galaxies with 
z < 0.14, My < —20.5 and Mbh > 10 6 Mq are shown. The blue points represent galaxies whose SMBHs accrete via an ADAF. These include galaxies 
experiencing both the radio and quasar modes. The red points represent galaxies powered by thin discs around SMBHs. The green points denote super- 
Eddington objects. The observational data are taken from Sikora et al. (2007): BLRGs are shown by open triangles; radio-loud quasars by filled triangles; 
Seyfert galaxies and LINERs by open squares; FRI radio galaxies by open circles; and PG Quasars by open stars. 



—2 regime. When compared to the data from Sikora et al. (2007), 
the predictions of the chaotic model give an adequate representation 
of the AGN on the 3? — A plane. The model reproduces well the 
loci of the different AGN populations and the slope of the overall 
distribution. Objects such as quasars and Seyferts are radiating at 
~ (0.01 — l)iEdd, while the radio galaxies are characterised by 
very low sub-Eddington luminosities (< 0.011/Edd)- 

By contrast, in the prolonged model there is an absence of ob- 
jects with log 3? < 3 in the —4 < A < 3 regime. This is mainly 
because most objects pile up in the upper envelope of the distribu- 
tion due to their high radio luminosities (and BH spin values). We 
note again that FR-I galaxies (open circles) have higher values of 
radio loudness than the models. This could be due to the fact that 
the measured radio luminosities for these objects include power 
from the lobes. 

The objects in our sample populating the radio-loud part of 
the JL — A plane are associated mainly with giant ellipticals that 
host very massive BHs (see below). Accretion onto these BHs in 
the ADAF regime produces significant radio power because the 
spin is very high. In combination with the low optical luminosi- 
ties of these objects, this gives rise to very high 31 values. There- 
fore, these specific properties of the host of the radio-loud AGN 
may suggest a correlation between 3J and Mbh- However, when we 
plot 3? versus Mbh (lower panels in Fig. I14t we find that in both 
models there is no apparent correlation between these two quanti- 
ties. The data display significant scatter along the Mbh axis which 
is due to the complex dependence of 3i on Mbh through the op- 
tical and jet luminosity expressions in Eqns. ( 144b . 045 1 . and l |35| >. 
Both models reproduce the loci of the different AGN population in 
the Sikora et al. data set remarkably well. ADAF systems, mainly 
identified with radio-loud sources, span a wide range of BH masses 
(10 6 — 10 10 Mq). In contrast, radio-quiet objects have BH masses 



in the range ~ 10 7 - 3 x 10 8 M©. Finally, in both models all 
the sources hosting BH with Mbh > 3 x 10 8 Mq are exclusively 
radio-loud. 



6.3 A physical view of the "chaotic accretion" population 

While both accretion models give a fair representation of the lo- 
cus of the heterogeneous radio-optical data points from Sikora et 
al. (2007), the chaotic accretion model provides a better fit to the 
much more homogeneous data from the radio luminosity function 
described in Section [54] Hence, in this section, we focus on the 
chaotic accretion model and explore the results in greater detail. 

In order to gain further insight into the predictions of the 
model, we show in Fig. Q3] how the main parameters that deter- 
mine the radio loudness of an AGN, a, Mbh and m, are distributed 
across the log Lb — log Lr plane for the chaotic accretion model. 

Comparing these plots to the data in Fig. |T3] the brightest 
radio-loud galaxies i.e. FR-I and BLRG, are identified with very 
massive (10 s — 10 9 M Q ), rapidly rotating BHs (spins > 0.8) ac- 
creting at the top of the ADAF branch at m ~ 0.01. In our galaxy 
formation model, these objects are hosted exclusively by giant el- 
lipticals and their central BHs grow through mergers and through 
the radio mode. The radio loud quasars (only 2 in our volume) are 
the extreme end of this population, and can be matched by the high- 
est mass objects in our sample. One issue with identifying these ob- 
jects with the ADAF branch is that they typically have high excita- 
tion spectra, showing that there is a bright UV disc (see e.g. March- 
esini, Celotti & Ferrarese 2004). However, the top of the ADAF 
branch is where the transition to a thin disc takes place. Observa- 
tions of stellar mass BH binary systems show that this transition 
is complex, probably taking on a composite structure with the thin 
disc replacing the hot flow at progressively smaller radii (see e.g. 
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Figure 14. Scatter plots of radio loudness (3J = (^b/^Sqhz) x ^r/^b) 
ra, BH mass (top panels) and A parameter (bottom panels) in the chaotic and 
prolonged accretion models for the sample in Fig. [13] (galaxies with z < 
0.14, My < -20.5 and M bh > 10 6 M o ). Different colours represent 
different accretion regimes: blue for ADAF systems and red for thin-disc 
systems. The observational data are taken from Sikora et al. (2007): BLRGs 
are shown by open triangles, radio-loud quasars by filled triangles, Seyfert 
galaxies and LINERs by open squares, FRI radio galaxies by open circles, 
and PG Quasars by open stars. 



the review by Done, Gieiiinski & Kubota 2007). Thus, better mod- 
elling of the details of the transition may well be able to reproduce 
the required UV ionising spectrum. 

The Palomar-Green (PG) Quasars are mostly concentrated 
at the top end of the narrow thin disc sequence, extending from 



10« to ■ 



10 45 erg s 1 in optical luminosity and from < 



to ~ 10 40 erg s _1 in radio luminosity. The higher mass accretion 
rates required to accrete via a thin disc are only generally possi- 
ble in lower mass BHs, 10 7 — 10 s M (see Fig.[4]l, so these have 
similar optical luminosities to the top of the ADAF branch as their 
higher mass accretion rate is partially cancelled by the lower mass 
BHs. These lower mass BHs are hosted in lower mass galaxies, 
predominantly spirals, and have lower spin, but the majority of the 
suppression in radio power comes from the ~ 2 orders of magni- 
tude stronger jet emission in the ADAF regime than in the thin disc 
regime at rh = 0.01. However, a few of the optically brightest radio 
quiet objects are powered by Eddington (or even super-Eddington) 
accretion onto 10 s Mg BHs, and these could be hosted by ellipti- 
cals. 

Seyfert and LINER galaxies constitute an interesting sample 
in our model because of their heterogeneity. Those exhibiting high 
nuclear luminosities (Lb > 10 43 ergs -1 ) (mainly Seyferts), ap- 
pear to occupy both the ADAF and thin-disc regimes, whereas 
those with lower luminosities (LINERs) are generally powered by 
an ADAF. These are the most numerous type of AGN seen in the 
local Universe. The majority of these objects have BHs with masses 
of 10 7 — 10 8 Mq. Thus, they have low-to-moderate spin and are 
relatively quiet at radio luminosities. These are the same population 
of objects as the PG Quasars, but at lower mass accretion rates. 
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Figure 15. Scatter plot of radio luminosity vs. B-band nuclear luminosity 
for the chaotic accretion model (as in Fig. l 13t . with model points coloured- 
coded according to BH spin (top), mass (middle) and accretion rate (bot- 
tom). The colours represent different ranges of values as indicated in each 
key. 
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7 DISCUSSION 

The process of BH growth in our model is dominated by the ac- 
cretion of large amounts of gas that allow a BH to double its mass 
several times during a Hubble time. The evolution of the spin of the 
BH is strongly influenced by the nature of the accretion episodes. 
It is often assumed that BHs accrete gas via a disc with constant 
angular momentum. In this case (the prolonged accretion model), 
BHs are systematically spun up during accretion, with most of them 
ending up with maximal spin. Even minor mergers or disc instabil- 
ities, typical of the growth of the bulge of spiral galaxies, can trig- 
ger a gas flow onto the BH of mass comparable to that of the BH 
itself and thus lead to a = 1. By contrast, if the accreting material 
fragments at its self-gravity radius and the associated star forma- 
tion randomises the angular momentum direction of each of the 
fragments (the chaotic accretion model), then accretion proceeds 
in multiple randomly oriented episodes. This results in a low spin 
BH. However, for the most massive BHs, the major growth chan- 
nel is not accretion but BH-BH mergers, and these generate fairly 
rapid spins of 0.7 — 0.8. Thus, in the prolonged accretion model, 
the most massive BHs have slightly lower spin than the bulk of the 
BH population, whereas in the chaotic model they have larger spin. 

Our results for the spin distributions are in good agreement 
with those of Berti & Volonteri (2008). However, our predictions 
for the prolonged-accretion model deviate significantly from the 
results of Lagos et al. (2009). These authors obtain a strong spin bi- 
modality when they consider accretion of gas via an accretion disc 
of constant angular momentum. The difference between the two 
approaches seems from the different amounts of cold gas available 
for accretion in the two models (see Kim et al. 2010 for a recent 
study of the cold-gas abundance in the Bower et al. 2006 model). 
In the original Bower et al. (2006) model and in the updated ver- 
sion that we have adopted in this study, the typical amount of cold 
gas accreted by a 10 6 — 10 s M© BH is enough to spin it up to 
the maximum value. This is not true, however, for the Lagos et 
al. model where cold gas accretion includes several episodes with 
M gEL3 < Mbh that do not spin up the BH efficiently. The difference 
in the final spin distributions illustrates the impact of the underly- 
ing galaxy formation model on the inferred properties of the BHs 
at the centre of galaxies. 

We have used the resulting distributions of BH mass, spin and 
mass accretion rate to calculate the properties of AGN, assuming 
that any accretion at a rate m < 0.01 proceeds via an ADAF, while 
accretion at higher rates forms a standard thin disc. The optical 
luminosity is then fairly straightforward to calculate, but the radio 
luminosity depends on the jet model. This is poorly understood, but 
both observations and theoretical models agree that the collapse of 
the thick ADAF into a thin disc leads to a similar collapse of the 
jet emission by several orders of magnitude. Thus, there is already 
a clear accretion mode switch in radio-loudness predicted by these 
models (Maccarone, Gallo & Fender 2003; Jester 2005). 

The main issue then is how jet power couples to spin. In the 
prolonged accretion model, there is a very small range in spin since 
all BHs are spinning rapidly. Thus, independently of the detailed jet 
model, the slope of the radio luminosity function is directly given 
by the relative numbers of low and high mass BHs accreting in 
the ADAF regime. This is steeper than observed since lower mass 
BH are much more numerous than higher mass ones. Instead, in 
the chaotic accretion model, the observed radio luminosity function 
can be reproduced if the jet power depends strongly on spin. The a 2 
dependence on spin power in the BZ jet models is sufficient to make 
a large difference in radio luminosity between the numerous low 



mass, low spin BHs and the much rarer high mass, high spin BHs. 
This model reproduces the shape of the observed radio luminosity 
function very well. 

A mass-spin correlation was also suggested by Sikora et 
al. (2007) to explain the range in the radio-to-accretion disc lu- 
minosity ratio seen in their (very heterogeneous) sample of AGN. 
However, the model they propose to establish this correlation is 
rather different from ours which is based on an actual calculation 
of galaxy formation from CDM initial conditions. Instead, Sikora et 
al. speculated that the high spin resulted from major mergers which 
triggered large gas flows with constant angular momentum direc- 
tion onto the nucleus, spinning up the BH and producing an ellipti- 
cal galaxy. Conversely, they argued that low spin resulted from mi- 
nor mergers of randomly aligned satellite galaxies (a random walk 
with spin up, spin down). These leave the galactic gas disc intact, 
producing a spiral galaxy. However, as discussed above, since most 
of the mass comes from the disc of the host galaxy, even minor 
mergers result in enough gas flowing to the centre to spin the BHs 
up to maximal. Fragmentation and chaotic accretion are necessary 
in order to produce low-spin BHs. Major mergers are indeed the 
key to the high spin of the most massive BHs, but this is the result 
of BH-BH mergers, not of gas accretion. 

One of the distinguishing features of our work is that the prop- 
erties of the BHs and their associated AGN are calculated ab initio 
within the context of a full model of galaxy formation in a ACDM 
universe. This model has been shown to agree with a large vari- 
ety of observational data such as galaxy luminosity functions in 
various passbands and at different epochs, galaxy colours, the cos- 
mic star formation history and structural scaling properties such 
as the Tully-Fisher relation and the Mbh — a relations (Bower et 
al. 2006). In this work we have augmented the galaxy formation 
model with calculations of BH spin and of the optical and radio 
output during accretion onto the BH. The two different accretion 
models we have explored (prolonged and chaotic) span the range 
of likely spin distributions and are good templates to investigate 
how the optical and radio luminosities of accreting systems depend 
on BH spin. As shown in Sections 5 and 6 and discussed in the 
previous paragraphs, the chaotic model gives a reasonable overall 
match to the observations, suggesting that a model in which BH 
spins develop a bimodal distribution is the most plausible. 



8 A UNIFICATION SCHEME FOR THE AGN ACTIVITY 

Finally, we present an interpretation of the optical and radio sig- 
natures of an active nucleus at different evolutionary stages of the 
host galaxy within a chaotic accretion framework. This is illustrated 
schematically in Fig. [16] where we show the relative position of 
the different types of active galaxies on the fundamental parame- 
ter (a, rh and Mbh) plane (left panel; the a and Mbh axes have 
been merged into one since the two quantities correlate as shown 
in the right-hand panel of Fig. |9]l and the jet/accretion-luminosity 
plane (right panel). We also speculate on the outcome of super- 
Eddington accretion. Such objects are rare in the local Universe, 
but are increasingly important at higher redshift (see Fig.[5J- 

We start by considering the sources with the lowest mass BHs. 
During a minor galaxy merger, cold gas from the disc is transferred 
to the bulge, triggering star formation and BH growth. This process 
represents a natural mechanism for the gradual growth of spheroids 
in spiral galaxies (Parry et al. 2009). Accretion of gas onto the BH 
usually occurs at sub-Eddington rates through a thin disc which 
turns the host spheroid into the bright nucleus of a Seyfert galaxy. 
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Figure 16. Left: the location of the different AGN types on the fundamental parameter (a, rh, M^h) plane. The a and Af^h axes have been merged into one 
since the two parameters are correlated as shown in the right hand plot of Fig. [9] in the chaotic accretion model. The shaded areas represent the different 
accretion regimes indicated by the keys. Right: the location of the different AGN types on the optical - radio plane. The shaded areas represent different 
accretion regimes. The arrows show how the physical parameters vary on the plane. Radio galaxies and LINERs lie in the ADAF regime. Seyferts are 
preferentially found in the thin-disc regime, even though a substantial population of Seyferts may be powered by an ADAF. RQQs accrete at relatively higher 
accretion rates than Seyferts and RLQs lie at the top of the ADAF branch. We speculate that some of the latter may be found also in the super-Eddington 
regime. 



The bright UV flux from this disc produces a strongly ionising 
spectrum, resulting in strong emission lines. However, gas flows 
close to the disc may also contribute to the emission (or absorp- 
tion) spectrum of the nucleus. The orientation of the central engine 
relative to a distant observer may have an important impact on the 
observed spectral features of the source and could account for the 
Type 1 and 2 sub-classification of Seyferts. 

If the host galaxy experiences a major merger or a galactic 
disc instability then the entire galactic gas disc is assumed to lose 
most of its angular momentum, participate in a starburst and is 
added to the stellar spheroid mass. This supplies the central region 
of the galaxy with large amounts of cold gas that feed the SMBH 
with several solar masses of material per year at a near- or super- 
Eddington rate. Thus, the nucleus becomes exceptionally luminous, 
resulting in a quasar. The collapse of the cold gas reservoir also ini- 
tiates an intense starburst which could be contemporaneous with 
the quasar phase. However, during the quasar phase, the galaxy is 
normally too dim to be seen against the vast amounts of radiation 
produced by the central engine. The radio output of the central en- 
gine when the accretion rate becomes super-Eddington is not clear 
since our models do not extend to that regime. Perhaps if the mag- 
netic field strength close to the BH is enhanced due to a transition 
to a thick disc when the flow exceeds the Eddington limit, strong 
jets may be launched establishing the galaxy as a RLQ. 

The duration of quasar activity is usually a few hundred mil- 
lion years. At the end of it (and the associated starburst), the system 
will have consumed or ejected its cold gas leaving a stellar bulge. 
In most halos, the ejected gas falls back into the galaxy and a new 
galactic disc is formed on a dynamical timescale. Further gas accre- 
tion, or galaxy mergers, may subsequently trigger another period of 



quasar activity and initiate a new growth era for the SMBH. Even- 
tually, the halo may grow massive enough that gas is no longer 
able to cool rapidly and the halo enters the hydrostatic cooling 
regime. In this case, the end-product may be an elliptical galaxy 
in which further gas cooling is restricted by the radio-mode feed- 
back. By construction, this accretion onto the BH takes place in the 
ADAF regime and the host elliptical appears as a radio galaxy. In 
the extreme case where the accretion occurs at the top of the ADAF 
branch and the central SMBH is rapidly rotating the galaxy is iden- 
tified as a RLQ. The orientation of the central engine and jets may 
explain the different subcategories of RLQs. For example, if the jet 
axis lies close to the line-of-sight of a distant observer, the source 
may be visible as a blazar. 



9 CONCLUSIONS 

We have presented a model of AGN activity in which the evolution 
of the galaxy is calculated in its full cosmological context and the 
luminosity of AGN at radio and optical wavelengths is calculated 
using a model for gas accretion and the generation of jets. We first 
considered the evolution of SMBH spin. We have found that the dif- 
ferent astrophysical processes that influence the growth of SMBHs 
have a significant effect on the global spin distribution. For exam- 
ple, if accretion of gas of constant angular momentum dominates 
the growth of SMBHs, the associated holes will be rapidly rotating. 
However, if the gas that is fed into the SMBH has random angular 
momentum, a bimodal spin distribution results. In this case, high 
spin values occur only for the most massive BHs (> 10 8 M©), 
and these are mainly due to gas poor major mergers, where the BH 
growth is dominated by the BH-BH merger. 
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We have coupled this mass, spin and mass accretion rate evo- 
lution to a model for the accretion flow and jet. The accretion flow 
is assumed to form a geometrically thick, hot, radiatively inefficient 
flow (ADAF) for rh < 0.01, and to collapse to a thin disc at higher 
mass accretion rates. The jet power couples strongly to the accre- 
tion mode since it depends on the vertical (poloidal) magnetic field 
component close to the BH horizon. The collapse by two orders of 
magnitude in the scale height of the flow results in a similar drop in 
radio power. This already produces a dichotomy in radio properties 
which explains the distinction between radio-loud and radio quiet 
objects. 

However, we also find that for our model to match the slope 
of the observed radio luminosity function, it is necessary for low 
mass BHs to generate relatively less jet power than high mass BHs. 
This is readily achieved in models where the jet couples strongly 
to spin, such as in the classic BZ jet mechanism, and where lower 
mass BHs have lower spin than the most massive BHs, as in our 
chaotic accretion model. 

Coupling the chaotic accretion model to the BZ jet mechanism 
results in an AGN population which reproduces the diversity of 
nuclear activity seen in the local Universe. In particular, the model 
accounts for the radio and optical luminosities of the FR-I, BLRG, 
Seyfert and LINER galaxy populations. This is the first consistent 
demonstration that a great part of the phenomenology of AGN can 
be naturally explained by the coeval evolution of galaxies and BHs, 
coupled by AGN feedback, in a cold dark matter universe. In future 
work, we will address the evolution of AGN across cosmic time, as 
a crucial test of galaxy formation models. 
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